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Abstract 

Thermoelectric  elements  made  of  semiconductor  plates  laminated  with  a  highly  electrically  and  thermally  conducting 
inter-connector  with  a  flow  channel  configuration  can  be  treated  as  integrated  Thermoelectric  Device  (TED). 
Elements  with  constructed  bulk  crystalline  n-type  and  p-type  (Bismuth-Telluride-Selenium)  semiconducting 
materials  and  copper  as  a  conducting  material  are  considered.  In  this  study,  the  thermoelectric  performance  of  such  an 
element  using  fluid-thermo-electric  coupled  field  numerical  methods  has  been  investigated.  The  TED  is  subjected  to 
constant  cold  temperature  at  the  bottom  and  top  surfaces  while  the  inter-connector  channel  walls  are  exposed  to  hot 
fluid  flow;  the  remaining  surfaces  are  kept  adiabatic.  The  performance  of  the  integrated  TED  element  is  studied  in 
terms  of  power  output  P0  ,  heat  input  Qh ,  conversion  efficiency  rj  and  produced  electric  voltage  V  for  different  load 
resistances  RL ,  inlet  hot  fluid  temperatures  Tin  and  flow  rates  ReDh  .  For  a  fixed  Tin ,  optimum  77  is  shown  at  a  load 
resistance  RL  which  is  slightly  lesser  than  total  internal  resistance  value.  An  increase  in  Tm  results  in  an 
enhancement  in  P0  and  77 ,  however,  it  has  a  minimal  effect  on  the  variation  of  optimum  load  resistance  RoptmL  .  At 
higher  Tin  values,  the  increment  in  RL  showed  a  significant  change  in  the  heat  input  Qh  values.  The  ReDh  had  a 
prominent  effect  on  TED  performance  predictions.  The  P0  increased  nearly  three-fold  at  Rem  =  500  in  comparison 
to  ReDh  =  100 .  It  is  recommended  that  for  an  accurate  modeling,  design  and  optimization  of  state-of-the-art  TED 
with  flow  channels  to  be  carried  out  using  fluid-thermo-electrical  coupling  field  simulations. 

©  201 1  Published  by  Elsevier  Ltd.  Selection  and/or  peer-review  under  responsibility  of  the  organizing 
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1.  Introduction 

Thermoelectric  device  (TED)  is  made  by  joining  two  dissimilar,  electrically  conductive  materials.  TED 
operates  based  on  the  principles  of  Seebeck  and  Peltier  effects.  Utilizing  the  Seebeck  effect,  TED  acts  as 
a  power  generator  when  the  two  junctions  are  exposed  to  temperature  differential.  The  TED  then 
generates  a  voltage  potential,  known  as  the  Seebeck  voltage.  Similarly,  TED  works  as  a  refrigerator  when 
electric  current  applied  across  the  element  creates  a  temperature  differential  at  the  material  junctions, 
known  as  the  Peltier  effect  [1],  These  devices  have  neither  moving  parts  nor  emit  noise  pollution  and  are 
also  environmental  friendly.  TEDs  have  a  wide  range  of  applications  include  temperature  measurement, 
power  sources  in  oil  and  natural  gas  facilities,  remote  radio  and  satellite  power  generation  stations, 
refrigeration  cooling  and  waste  heat  recovery  systems. 

It  is  estimated  that  nearly  70  %  of  the  fuel  energy  in  automobiles,  power  plants  and  industrial  systems  is 
rejected  to  surroundings  as  a  waste  heat.  A  reliable  technology  for  effectively  converting  waste  heat  into 
useful  energy,  i.e.  electricity  is  TEDs.  Conventional  cascading  and  segmented  TEDs  are  the  primary 
methods  for  extracting  waste  heat  from  low  and  high  temperature  sources.  In  [2],  they  developed  a 
segmented  TED  using  novel  p-type  Zn4Sb3  and  CeFe4Sb14,  and  n-type  CoSb,  based  alloys,  which 
operates  over  a  temperature  range  of  300-973  K.  They  achieved  a  conversion  efficiency  of  15  %  . 
Recently,  [3]  studied  using  a  one-dimensional  thermal  resistance  model  for  thermoelectric  modules  with 
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applications  to  waste  heat  recovery  from  an  automobile  engine  and  compared  their  model  with 
experimental  data.  The  results  showed  that  the  performance  of  a  module  on  the  exhaust  pipe  performs 
better  compared  to  a  module  on  the  radiator  system. 

Most  of  the  studies  dealt  with  one-dimensional  models,  even  though  electric  field  and  temperature 
distributions  are  multidimensional  in  TEDs.  However,  in  a  few  articles  [4,5],  they  demonstrated  that  the 
three-dimensional  simulations  are  more  accurate  in  modeling,  design  and  optimization  of  TEDs.  The 
authors  [4]  showed  that  the  electric  power  produced  is  sensitive  to  hot  fluid  properties  and  significantly 
depends  on  real  operating  conditions,  but  at  present  is  inefficient  if  one  considers  the  weight  to  power 
ratio.  In  [5],  they  proposed  and  implemented  a  non-linear  fluid-thermal-electric  three-dimensional 
numerical  model  for  thermoelectric  generators  in  FLUENT-UDS  (User  Defined  Scalar)  environment. 
Although  conventional  TED  designs  are  applicable  to  waste  heat  recovery  operations,  they  require  large 
amount  of  rare  earth  materials  for  module  fabrication.  These  designs  also  induce  large  amounts  of  thermal 
stresses  due  to  their  structure.  Furthermore,  there  exist  substantial  thermal  resistances  among  the 
electrically  conductive  copper  shunts,  the  typical  ceramic  plates  and  heat  sinks.  Therefore,  keeping  in 
mind  the  above  salient  points,  a  novel  integrated  TED  structure  is  proposed,  where  a  major  part  of  the 
semiconductor  material  is  replaced  with  a  conductor,  such  as  copper,  along  with  flow  channels  drilled 
through  it.  This  flow  channel  configuration  acts  as  a  low  thermal  resistance  heat  exchanger  eliminating 
the  thermal  resistances  attributed  to  the  heat  sinks.  This  novel  design  helps  reduce  thermal  stresses, 
materials  usage  and  the  thermal  resistances  without  compromising  its  conversion  efficiencies  and 
performance.  In  this  paper,  a  fluid-thermo-electrically  coupled  three-dimensional  numerical  model  for 
TEDs  is  implemented  in  the  FLUENT-UDS  package.  Furthermore,  using  developed  numerical  model,  the 
performance  of  such  a  novel  integrated  TED  is  investigated  for  different  load  circuitry  resistance,  hot 
fluid  inlet  temperature  and  flow  rate  values. 

2.  Geometry,  Governing  equations  and  Boundary  conditions 

The  schematic  of  a  novel  integrated  Thermoelectric  Device  (TED)  with  flow  channel  configuration  being 
investigated  is  shown  in  Fig.  1 .  The  integrated  TED  element  with  square  cross-section  WxD  consists  of 
a  copper  inter-connector  laminated  between  crystalline  n-type  (75%  Bi2Te:  25%  Bi2Se3 )  and  p-type 
(25%  Bi2Te3  75  %  Sb2Te2,  1.75%  excess  Se  )  bulk  materials.  The  inter-connector  rectangular  flow 
channel  size  is  0.8(7,  -  2d)  x  0.6D  .  The  top  and  bottom  surfaces  of  the  TED  element  are  connected  to 
copper  connectors.  The  end  surfaces  of  these  connectors  act  as  electrical  terminals  for  load  circuitry  of 
load  resistance  R,  .  The  top  surface  of  the  upper  connector  and  bottom  surface  of  the  lower  connector  are 
subjected  to  constant  cold  temperature  Tc .  A  hot  fluid  with  constant  temperature  Tin  and  uniform  velocity 
Uin  enters  the  main  flow  channel  of  cross-sectional  size  (Z,  -  2  d)xD  at  the  inlet,  flows  the  through  inter¬ 
connector  channel  and  leaves  at  the  outlet  (as  shown  in  Fig.  1).  The  inter-connector  left  and  right  side 
surfaces  and  it’s  channel  walls  are  exposed  to  the  hot  fluid.  The  remaining  surfaces  of  TED  and  the  main 
flow  channel  walls  exposed  to  ambient  are  subjected  to  adiabatic  condition. 


2090 


B.  V.  K.  Reddy  et  al.\/ Energy  Procedia  14  (2012)  2088  -  2095 


Fig.  1.  Schematic  of  integrated  thermoelectric  device 


The  fluid  flow  and  heat  transport  in  the  fluid  domain  are  governed  by  the  continuity  of  mass,  momentum 
and  energy  equations  subjected  to  certain  assumptions.  The  flow  is  steady,  laminar  and  incompressible. 
Thermo-physical  properties  of  the  fluid  are  invariant  with  temperature.  In  the  thermoelectric  material,  the 
current  flow  and  energy  transport  are  governed  by  the  continuity  of  current  density  and  the  energy 
equations  subjected  to  following  assumptions:  the  materials  are  heterogeneous  and  isotropic  with 
temperature-dependent  thermoelectric  properties;  electrical  resistivity  p,  Seebeck  coefficient  a  and 
thermal  conductivity  k,  are  specified  as  polynomial  functions  of  temperature  [6];  the  thermal  and 
electrical  contact  resistances  at  the  interface  of  n-type  or  p-type  and  copper  materials  are  negligible.  The 
system  of  partial  differential  equations  governing  the  fluid-thermo-electric  phenomenon  in  TEDs  under 
steady-state  conditions  for  the  continuity,  momentum  and  heat  transport  [1]  is  written  as 
For  fluid  domain:  V.v  =  0  (1) 

py(v.Vv)  =  -VP+  //V2v  (2) 


pfcp  (v.  V  T)  =  V.(kfV  T) 


(3) 


For  thermoelectric  material: 


V.J  =  0 

V(*v:r)  +  pi2  -  TJ. ^ j' VT  +  (Va)T  J  =  0 


(4) 

(5) 


where  J  is  the  current  density  and  it  is  evaluated 


J 


P 


[VV  +  aVr]=  0 


(6) 


In  Eq.  (5),  the  second  and  third  terms  on  the  left-hand  side  represent  the  Joule  heating,  Peltier  and 
Thomson  effects  respectively.  In  Eq.  (6),  the  first  and  second  terms  on  right-hand  side  are  the  current 
densities  due  to  electrostatic  potential  distribution,  and  the  Seebeck  electric  potential  in  the  thermoelectric 
material,  respectively. 
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The  associated  boundary  conditions  for  Eqs.  (1)  to  (6)  with  respect  to  the  geometry  shown  in  Fig.  1  are: 

At  the  terminal  in:  J  =  —  = - — -  and  —  =  0  (7) 

4  4(4  +  4) 

where  V0  is  the  total  build-in  open  circuit  voltage  generated  in  the  element  at  a  no-load  condition,  which 
is  known  as  the  Seebeck  voltage.  Rln  is  the  total  internal  resistance  equaling  to  the  summation  of  n-type, 
p-type,  connectors  and  inter-connector  materials  resistances  and  R,  represents  the  external  load 
resistance. 


At  the  terminal  out:  V  =  0  and  —  =  0 

n 

At  the  top  wall  of  the  upper  connector  and  bottom  wall  of  the  lower  connector: 


(8) 


(9) 


All  other  surfaces  of  the  TED  are  exposed  to  surroundings:  —  =  0 

At  the  fluid  flow  inlet:  x  =  0 :  u  =  Um,v=w=  0  and  T- 

At  the  fluid  flow  outlet:  x  =  (L\  +  W  +  L2) :  -^-4  =  =  v  =  w  = 

d df 


0T 

5# 


=  0 


(10) 

(11) 

(12) 

(13) 


At  the  channel  walls  exposed  to  surroundings:  u  =  v  =  w=  0  and 

where  <+  denotes  the  direction  normal  to  the  corresponding  boundary  wall.  L\  =  5W  and  L2  =  20W  are 
the  lengths  of  the  upstream  and  downstream  boundaries  of  the  TED,  respectively  (as  shown  in  Fig.  1).  At 
the  interface  between  the  n-type  or  p-type  and  conductor  materials,  the  continuity  of  current  density, 
temperature  and  heat  flux  conditions  are  imposed. 


The  power  output  from  the  TED  for  a  given  load  resistance  RL  ,  the  total  heat  input  at  the  inter-connector 
walls  exposed  to  flowing  fluid  and  the  thermal  conversion  efficiency  are  evaluated  as: 


P0  =  I2Rl 


and  Q.  =  -kf  f  —  dA  rj  =  — 

h  s  a 


(14) 


3.  Numerical  solution  methodology  and  validation 

The  numerical  simulations  are  done  using  the  finite  volume  formulation  of  Eqs.  (1)  to  (5)  and  the 
associated  constitutive  relation  (6)  along  with  boundary  conditions  Eqs.  (7)  to  (13)  in  ANSYS  FLUENT- 
UDS  (User  Defined  Scalar)  environment.  The  integrated  TED  model  and  mesh  are  generated  in  Gambit 
2.4.  The  pressure-velocity  coupling  is  handled  using  SIMPLE  algorithm.  The  convective  term  in  Eq.  (2), 
diffusion  terms  in  Eqs.  (2)  and  (4)  and  pressure  term  in  Eq.  (2)  are  descritized  with  a  power  law  scheme. 
The  Seebeck  potential  and  Ohmic  electric  potential  (Eq.  (6))  are  calculated  using  UDS  fields.  The  Ohmic 
heating,  Peltier  and  Thomson  effects  are  modeled  as  source  terms  in  the  energy  equation  (Eq.  (5)).  The 
electric  current  is  evaluated  based  on  the  open-circuit  Seebeck  voltage  ( V0 )  produced  at  the  given  load 
resistance  RL  (Eq.  (7)).  The  detailed  numerical  implementation  has  been  demonstrated  in  [5], 


Table  1 .  Comparison  of  thermoelectric  generator  parameters  Th=42TC  and  Tc=2TC 
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Quantity 

Analytical  [61 

|  ANSYS  (a)  [71 

|  ANSYS  (b)  [71 

|  Present  (a) 

|  Present  (b) 

Qh,W 

13.04 

13.03 

11.07 

13.0 

10.72 

P0,W 

1.44 

1.43 

1.05 

1.43 

1.07 

q,  % 

11 

11 

9.5 

11 

9.98 

I,  A 

19.2 

19.1 

16.4 

19.09 

16.5 

The  convergence  criteria  (the  difference  between  consecutive  iterations  over  a  domain)  set  for  the  flow  is 
10”4  and  for  the  current  density,  energy,  Seebeck  and  Ohmic  electric  potentials  are  10"12  .  Grid 
independence  tests  are  performed  for  integrated  TED  (as  shown  in  Fig.  1).  A  grid  size  622800  of  non- 
uniform,  orthogonal  grid  used  in  the  simulations.  The  developed  three-dimensional  thermoelectric 
numerical  code  has  been  validated  with  published  results  of  a  thermoelectric  generator  for  the  cases  of  (a) 
constant  properties  and  (b)  temperature-dependent  properties  as  given  in  Table  1.  The  results  are  in  good 
agreement  with  the  published  both  analytical  and  numerical  results. 

4.  Results  and  Discussion 

The  numerical  simulations  are  conducted  on  a  integrated  thermoelectric  device  (TED)  to  study  its 
performance  when  coupled  with  fluid  flow  phenomenon.  The  geometric  dimensions  i.e.  cross-sectional 
area  (W  =  5  mm  and  D  =  5  mm)  and  height  (d  =  5  mm)  of  the  n  or  p-type  element,  the  connector 
(2W x  Dx  1.5),  the  interconnector  ((L-2d)  x  D)  and  its  channel  (0.3Dx  0.8(L-2d))  sizes  are  invariant.  For 
fixed  cold  wall  temperature  Tc  =  300  K,  the  performance  of  the  TED  in  terms  of  the  power  output  P0 , 
heat  input  Qh,  conversion  efficiency  q  and  Ohmic  electric  voltage  V  for  various  hot  fluid  inlet 
temperature  350  <  Tin(K)  <  550,  the  load  circuitry  resistance  10-4  <  RfOhms)  <  HE1  and  the  Reynolds 
number  based  on  the  hydraulic  diameter  Dh  50  <  Rem  <  500  has  been  investigated  in  detail. 

4.1.  Effects  of  load  resistance  RL 

For  fixed  Rem  =  500  ,  Fig.  2  shows  the  variation  of  power  output  P0,  heat  input  Qh,  and  conversion 
efficiency  q  and  Ohmic  voltage  drop  V  with  load  resistance  R,  for  different  inlet  hot  fluid 
temperatures  Tjn .  From  Figs.  2a  and  2c,  for  a  given  Tin ,  there  exists  a  maximum  power  output  P0max  when 
Rl  is  equal  to  Rw  and  the  TED  has  optimum  thermal  efficiency  q  at  R,  =  RoptmL ,  typically,  which  is  not 
equal  to  Rm  .  The  optimum  load  resistance  RoptmL  increases  marginally  with  an  increase  in  Tm  values 
while  all  other  are  parameters  kept  constant.  It  is  noticed  that  the  Rin  value  is  always  slightly  lesser  than 
Rop,m,  for  a  given  Tin  condition.  At  higher  Tin  values,  the  increment  in  R,  showed  significant  variations 
in  Qh  predictions,  however,  a  minimal  decrease  at  lower  values  of  Tin  is  observed  (in  Fig.  2b). 

From  Fig.  2d,  for  a  given  Tin  value,  an  increase  in  R,  results  in  a  decrease  in  V  values.  This  is  due  to  the 
gain  in  total  resistance  Rin  +  RL  value  of  the  TED.  On  the  other  hand,  at  a  given  R,  value,  an  increment 
in  Tjn  shows  an  enhancement  in  V  values.  This  is  due  to  raise  in  the  Seebeck  voltage  through  the 
temperature  differential  between  the  inter-connector  channel  wall  and  the  cold  surfaces  ( Th  -  T )  . 
Furthermore,  irrespective  of  Tin  value,  the  responses  of  V  with  RL  show  same  behaviour. 
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Fig.  2  The  performance  parameters  variation  with  load  resistance  ( RL  )  (a)  Electrical  power  output  (b)  Heat  input  (c)  Thermal 
efficiency  and  (d)  Ohmic  voltage  for  various  TjK  values  [at  ReDh  =  500,  d  =  5  mm,  L  =  20  mm,  Tc  =  300  K] 


4.2.  Effects  of  Reynolds  number  Rem 


Figures  3  and  4  show  the  performance  of  a  integrated  TED  with  ReDh  values  for  different  Tm  values  at 
optimum  load  resistance  RoptmL  conditions.  For  a  fixed  Tin ,  as  the  inlet  flow  rate  Rem  increases,  an 
enhancement  in  the  performance  parameters,  namely  P0,  Qh,  q  and  V ,  is  observed.  It  indicates  that  the 
higher  Rem  improves  the  heat  transfer  rate  from  inter-connector  channel  walls  to  TED,  and  hence,  more 
temperature  differential  occurs  between  hot  and  cold  end  surfaces,  which  in  turn  produce  a  higher 
Seebeck  voltage  across  the  semiconductor  elements.  It  is  noticed  that  as  Tin  raises,  the  change  in  q  with 
Rem  becomes  prominent  and  the  behaviour  also  reflects  in  P0,  Qh,  and  V  values.  The  reason  could  be 
the  formation  of  recirculating  zones  near  the  inter-connector  walls  upstream  and  downstream  of  TED. 


In  Fig.  4a  and  b,  represent  the  x-velocity  and  temperature  contours  with  Tin  =  550 K  at  the  mid-plane  (as 
shown  in  Fig.  la)  for  Rem  =  50,  200  and  500  values,  respectively.  For  clarity,  a  downstream  length  of 
4 IT  and  a  upstream  length  of  2.5W  are  only  depicted  in  these  figures.  From  the  velocity  contours,  in  the 
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Fig.  3  The  variation  of  (a)  Electrical  power  output  (b)  Heat  input  (c)  Thermal  efficiency  and  (d)  Ohmic  voltage  with  Reynolds 
number  (  ReDh  for  different  Tin  values  at  Roplm  values  [at  d  =  5  mm,  L  =  20  mm,  Tc  =  300  K] 


Fig.  4  The  distribution  of  (a)  x-velocity  and  (b)  temperature  on  the  depth-mid-plane  at  different  ReDh  values  (i)  50  (ii)  200  and  (iii) 
500  [at  7,„  -  550  K,  d  -  5  mm.  I.-  20  mm,  Tc  =  300  K and  Ropm  =  6.5  x  1  O'3] 


downstream,  it  is  seen  that  the  flow  accelerating  towards  the  channel  walls  increases  as  the  inlet 
flow  Rem  raises.  This  flow  behaviour  helps  in  the  enhancement  of  heat  transfer  due  to  the  formation  of 
the  recirculation  zones  near  inter-connector  side  walls.  This  trend  has  been  revealed  through  raise  in  the 
inter-connector  channel  walls  temperature  as  shown  in  Fig.  4b.  The  larger  temperature  gradients  in  the 
thermoelectric  materials  facilitate  higher  P0  values  (as  shown  in  Fig.  3a). 


5.  Conclusions 

The  fluid-thermo-electric  coupling  mechanisms  such  as  the  Joule  heating,  Seebeck  and  Peltier  effects, 
Thomson  cooling/heating  and  convective  heat  transfer  in  thermoelectric  devices(TEDs),  along  with 
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temperature-dependent  material  properties  has  been  implemented  in  a  FLUENT-UDS  environment.  By 
focusing  on  improvement  in  power  output,  minimization  of  thermal  stresses  and  resistances,  and  reducing 
the  usage  of  rare-earth  materials  in  TED,  without  compromising  its  efficiency,  a  novel  integrated  TED 
has  been  proposed.  The  thermoelectric  performance  of  such  a  device  applied  to  waste  heat  recovery  in 
terms  of  power  out  P0 ,  heat  input  Q, ,  conversion  efficiency^  and  Ohmic  voltage  V  has  been  investigated 
using  numerical  simulations. 

For  a  fixed  inlet  flow  rate  Rem  and  fluid  temperature  Tm  ,  the  change  in  load  resistance  R,  results  in  an 
optimum  conversion  efficiency  q  and  maximum  power  output  P0max  at  R,  equaling  to  the  total  internal 
resistance  Rm  and  R,  equaling  to  the  optimum  load  resistance  R„ptm, ,  respectively.  The  R„ptm,  increases 
marginally  with  an  increase  in  Tin  values  is  noticed  when  all  other  parameters  kept  constant. 
Simultaneously,  a  decrease  in  V  values  with  raise  in  R,  values  is  observed. 

For  a  given  Rt  value,  an  increment  in  Tin  shows  an  enhancement  in  TED  performance.  It  is  also 
observed  that  at  higher  Tin  values,  the  increment  in  R,  showed  a  significant  change  in  Qh  predictions, 
however,  a  marginal  decrease  at  lower  values  of  Tjn . 

The  higher  flow  rate  Rem  enhances  the  heat  transfer  from  the  inter-connector  channel  walls  of  the  TED, 
and  hence  produces  a  higher  P0 .  The  formation  of  recirculation  zones  in  the  upstream  and  downstream  of 
TED  assist  in  heat  transfer  Qh  enhancement  with  an  increase  in  Tin .  Therefore,  the  fluid-thermo-electrical 
coupling  phenomenon  should  be  considered  for  detailed  analysis  of  TED  performance. 
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